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Abstract 

Based on a semiclassical expansion for quantum chromodynamics in the instanton liquid back- 
ground, the correlation function of the ++ scalar glueball current is given, and the properties 
of the ++ scalar glueball are studied in the framework of Gaussian sum rules. Besides the pure 
classical and quantum contributions, the contributions arising from the interactions between the 
classical instanton fields and quantum gluons are come into play. Instead of the usual zero-width 
y-^- \ approximation for the resonance, the Breit-Wigner form for the spectral function of the finite- width 

resonance is adopted. The family of the Gaussian sum rules for the scalar glueball in quantum 
chromodynamics with and without light quarks is studied. A consistency between the subtracted 
and unsubtracted sum rules is very well justified, and the values of the decay width and the cou- 
pling to the corresponding current for the ++ resonance, in which the scalar glueball fraction is 
M-f dominant, are obtained. 

D . PACS numbers: 11.15.Tk, 12.38.Lg, 11.55.Hx, 11.15.Kc, 12.39.Mk 

(N 
> 

in 

\o 
^t 

(N 

l> 

O 

o 



X 



To whom correspondence should be addressed. Email: jpliu@whu.edu.cn 



I. INTRODUCTION 



Glueballs, being composed of pure gluons in the limit of no quark fields, have at- 
tracted much attention since the theory of the strong interactions, quantum chromody- 
namics (QCD), was founded in the late of 1970s [IrM]- Estimates of glueball properties were 
obtained in a variety of approaches [5], ranging from model analyses |6l-ll3| , quenched lattice 
QCD (QLQCD)HJ|I3 and unquenched lattice QCD (UQLQCD) simulations 0, Q to 



QCD sum rule calculations [21H33 . 

The lowest ++ scalar glueball state is, in fact, the most intricate hadron state which 
is difficult to figure out. In lattice QCD (LQCD) the mass scale of the scalar glueball is 
predicted to be in the rang of 1.3-1.7 GeV |l4J-|l8J]. The ++ scalar resonances closest to this 



energy range are /b(1500) and /o(1710Hn the present data [341 ]. and some authors favor the 



former as the lightest scalar glueball [35|, ]36| , while some others prefer the latter [lOj, Il8 



Furthermore, both / (1500) and / (1710) may not be pure glueballs; to the contrary, these 



resonances can be considered to be the mixture of glueball and mesons [7H9J, |37J. In the 
QCD sum rule approach, the results of the mass of the scalar glueball are also different 
from each other. In the early days of QCD, some found a light ++ scalar glueball in the 
region of 300-700 MeV [3|, |30h33I| using the subtracted sum rule (SSR), while the others 
obtained a much heavier one in the 1-2 GeV region 38|, ]39j by using the unsubtracted sum 
rule (USSR). The inconsistency between both subtracted and unsubtracted sum rules has 
bothered scientists for many years. 

It should be noticed that, in the early QCD sum rule approach, it was already recognized 
that there is an onset of the nonperturbative physics (a departure from the asymptotic 
freedom) in the scalar glueball correlator at unusually short distances \x\ <£. A^py-. J4J. Such 

hard nonperturbative effects are usually considered to be coming from direct instantons, i.e. 
the tunneling processes which rearrange the QCD vacuum topology in localized regions 40, 



111. This physics was ignored in the early glueball sum rules except in Ref. [42]. 

Recently, a sizeable instanton contribution to the QCD sum rules of the ++ scalar 
lueball has been found on the basis of the instanton liquid model of the QCD vacuum 
4ll |. and supported by the lattice simulations [42j-]46|. However, the instanton-induced 



continuum contributions are neglected in the early QCD sum rule approach, and thus a 
reliable estimate of the glueball properties cannot be obtained. On the basis of the instanton- 
improved operator-product expansion (OPE), the authors in Refs. 12 31-1251 included the 



instanton-induced contributions to the continuum spectrum, made a great improvement 
to the consistency between different types of ++ scalar glueball sum rules, and gave new 
predictions for the mass and decay constant confirmed later by the Gaussian sum rule 



(GSR) calculation 27]. Moreover, we have clarified that the stability and the consistency 
for the SSR and USSRs for the ++ scalar glueball can be obtained by using a systematical 
semiclassical expansion of the background instanton fields in the instanton liquid model of 
the QCD vacuum \^. 

Instantons, which make a great difference in the scalar and pseudoscalar channels J4J, |4£ 
|49| . play a major role in the gluonic interactions, especially in the nonperturbative region 
J40|. Direct instanton contributions are included in the QCD sum rules J23l-l27j|. The results 
are found to be in good agreement with different USSRs. The compatibility between the 
SSR and USSRs become much better but are not satisfying yet, because only the leading 
classic effects are included in most of the calculations with the perturbative contribution 
and condensate contributions. It should be noted that the direct-instanton approximation is 



criticized due to the problem of double counting [241 ]. because both condensate and instanton 
contributions are included in the correlator, but most of the gluon condensates could be 
produced from instantons. 

Up to now, most of the theoretical calculations for the ++ scalar glueball were based on 
QCD Laplace sum rules (LSR), which emphasize the contributions of the lowest resonance, 
and have shown the power in the investigation of the nonperturbative properties of the 
hadron ground states. On the other hand, it should be noticed that the GSR developed later 
emphasizes only the contribution of the state considered, and has a cleaner background in 
comparison with the LSR. As a cross-check, the same problem should be investigated using 
GSR, because both Laplace and Gaussian sum rules are derived from the same underlying 
dynamical theory, and should give almost the same results. 

Motivated by the above considerations, our main objective in this paper is to investi- 
gate the ++ scalar glueball in the frame work of Gaussian sum rules. For the correlation 
function, we include the contributions from the interactions between the quantum gluons 
and the classical instanton background besides the ones comeing separately from these two 
different field configurations. For the spectral function, instead of using the usual zero- width 
approximation, we adopt the Breit-Wigner form for the considered resonance with correct 
threshold behavior, in order to get the information of not only the mass scale but also the 
full decay width. Moreover, without using the scheme of the mixture of the traditional 
condensates and the so-called direct-instanton contribution, we are working in the frame- 
work of the semiclassical expansion of QCD in the instanton liquid vacuum, a well-defined 
self-consistent procedure for the quantum theory justified by the path-integral quantization 
formalism. The paper is organized as follows: In Sec. [TT1 we give the expression for the 
correlation function of the ++ scalar gluon current. The spectral function corresponding 
to this current is constructed in Sec. IIHI Then, a family of the finite-width Gaussian sum 
rules are derived in Sec. IIVI In Sec. |V] the numerical simulation is carried out. Finally, a 
summery of our conclusions and a discussion are given in Sec. |VT1 



II. CORRELATION FUNCTION 

The correlation function for the scalar glueball in the Euclidean space-time with a virtu- 
ality q 2 is defined by 

n(g 2 ) = [d 4 xe l »- x (n\TO s (x)O s (0)\ty, (1) 

where \Q) is the physical vacuum, and O s the scalar glueball current with the quantum 
numbers J PC = ++ 

O s = a s G; u (B)G a ^(B), (2) 

in which, a s denotes the strong coupling constant. The scalar glueball current O s is gauge- 
invariant, and renormalization-invariant at one- loop level. In the spirit of the semiclassical 
expansion, and in order to maintain the 0(4) covariance, the gluon field strength tensor 
G^ U (B) is considered as a functional of the full gluon potential, B^ a = A^ a + a Ma , with A^ a 
and a^a being the instanton fields and the corresponding quantum fluctuations. 

The theoretical expression, n^ CD , for the correlation function II may be divided into the 
following three parts 

n QCD (Q 2 ) = n pcrt (Q 2 ) + n inst (Q 2 ) + n int (Q 2 ), (3) 



where Q 2 = q 2 , and Tl peTt (Q 2 ), IT nst (Q 2 ), and IT mt (Q 2 ) are the contributions from the only 
perturbative QCD, the pure instanton dynamics, and the interactions between the instantons 
and the quantum gluon fields, respectively. 

The perturbative contribution IP eTt (Q 2 ) is already known to be 



IP ert (Q 2 ) = Q 4 ln 
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where p? is the renormalization scale in the MS-dimensional regularization scheme, and the 
coefficients with the inclusion of the threshold effects are 
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for quarkless QCD up to two-loop level [5l|. Both expressions for n pert (Q 2 ) with and without 
quark loop corrections are used in our calculation for comparison. With the assumption that 
the dominant contribution to U mst (Q 2 ) comes from a Belavin-Polyakov-Schwartz-Tyupkin 
single instanton and anti-instanton solutions |52l-l54| and the multi-instanton effects are 
negligible (see a QCD spectral sum rule (QSSR) approach [251 ]) , and in view of the gauge- 
invariance of the correlation function, one may choose to work in the regular gauge of the 
classical single instanton potential 
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where i] a is the 't Hooft symbol, and Xq and p denote the position and size of the instanton, 
respectively. The pure instanton contribution is obtained to be [3|, |23|, [24J, |42|, [55 



n inst (g 2 
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where K 2 {x) is the McDonald function, n = J dpn(p) and p are the overall instanton 
density and the average instanton size in the random instanton background, respectively. 

It is noticed that the contribution to n^ from the interactions between instantons and 
the quantum gluon fields is of the order of the product of a s and the overall instanton density 
n. There is no reason to get rid of this contribution in comparison with the perturbative 
contributions of the higher order a^ considered in Il pert . To calculate such contribution, our 
key observation is that the instanton potential A a obeys also the fixed-point gauge condition 



(a; - x )„A a Jx - x ) = 



(9) 



due to the antisymmetry of the 't Hooft symbols. As a consequence, the instanton potential 
can be expressed in terms of the corresponding field strength tensor as follows: 

A a (x-x Q )= I duuF^[u(x - x )](x - x Q ) v (10) 

Jo 

and the gauge link with respect to the instanton fields is just the unit operator, and thus the 
trace of any product of the gauge-covariant instanton field strengths at different points is 
gauge- invariant. This allows us to conclude that the remainder quantum corrections to the 
gauge-invariant correlation function, arising from the interactions between the instantons 
and the quantum gluons, is gauge-invariant as well. Therefore, one may choose any specific 
gauge in evaluating the quantum correction to Il mt . 

It is noticed that the interference between the pure quantum part of the scalar glueball 
current, (O s ) quan tu m , and the pure classical instanton part, (O s )i nsta nton is vanishing due to 
the fact that the momentum integration of the massless gluon propagator or the product 
of massless quark and gluon propagators has no scale parameters, and is exactly zero in 
the dimensional regularization scheme. The first nonvanishing contribution comes from the 
contraction between the two quantum glueball currents with instanton legs. Working in 
Feynman gauge, our result for n mt is 

n int (Q 2 ) = C a s nir + a 2 s n[C x + C 2 (Qp) 2 \n(Qp) 2 + C 3 \n(Qp) 2 + C 4 (Qp)-% (11) 

where the coefficients are 

C = 62.62, C x = 1533.15, C 2 = 825.81, C 3 = -496.33, C 4 = -348.89. (12) 

It is remarkable to note that the fixed point Xq , which characterizes the gauge condition (TTU|) , 
disappears in the expression of Il mt , as expected from the gauge invariance of our procedure. 
The detail calculation for Il mt is much involved, and will appear elsewhere. 

III. SPECTRAL FUNCTION 

Now, we turn to specify the spectral function for the correlation function of the scalar 
glueball current. The imaginary part of the correlation function Eq. ([3]) is 



-TVS 2 



Imn QCD ( S ) = 

-16tc 4 s 2 np i J 2 (p v G)Y 2 (p v G) + a 2 s nn(C 2 p 2 s - C 3 ). (13) 
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The usual lowest resonance plus a continuum model is used to saturate the phenomenological 
spectral function, 

Imn PHE (s) = p had (s) + 6(s - s )Imn QCD (s), (14) 

where So is the QCD-hadron duality threshold, and p had (s) the spectral function for the 
lowest scalar glueball state. Instead of using the zero-width approximation as usual, the 
Breit-Wigner form for a single resonance is adopted for p had (s) in the quarkless world 



had,. PmT 



p*( s ) ( s - m 2 + r 2 /4) 2 + m 2 r 2 ' ^ 



where f 3 = (Q\O s \0 ++ ) is the coupling of the lowest resonance to the scalar glueball current 
Eq. ([2]). Recall the threshold behavior for p had (s) 



f 3 ->• A s, for s ->■ 0, 



(16) 



which is deduced from a low-energy theorem [4J, |38|, |56j . The early QCD sum rule approach 
had often used f 3 — > \qs in the whole lowest resonance region, however the obtained mass 
scale is too low to be expected in comparison with the lattice QCD results. In fact, the 
threshold behavior (TT6|) is only proven to be valid in the chiral limit; it may not be extrapo- 
lated far away. Therefore, instead of considering the coupling / as a constant [25fl , we choose 
a model for / as 

X Q s, for s < m^. 

X ml + A 3 , for s > ml ' 



f 



(17) 



where Ao and A are some constants determined late in numerical simulation, so that the 
spectral function p had (s) has the almost complete Breit-Wigner form with correct thresh- 
old behavior, and in cooperation with the threshold behavior which is important for the 
convergence of the corresponding integral. In Eq. (1T7|) . the constant A is invoked for the 
discontinuity at the chiral symmetry breaking. 

Although glueballs are well defined in quarkless QCD, the mixing with mesons makes the 
phenomenological side more complicated in full QCD; all scalar hadron states having a glue 
content should be interpolated by the gluonic current, and then present in the correlator 
with different couplings. In Refs. [7H9l. I34J-|36| . the glueball is shared between the three scalar 
hadrons. These three hadrons are very close in mass; the assumption of single resonance 
maybe be not appropriate. The form of the spectral function for three resonances we will 
use is taken to be 
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where all fi have the form shown in (|T7j) with the same Ao, because Ao is fixed to be 5 
GeV by the low-energy theorem of QCD, and so that its value is independent of what an 
individual resonance considered. 



IV. FINITE- WIDTH GAUSSIAN SUM RULES 



A family of Gaussian sum rules can also be constructed from the Borel transformation of 
the correlation function in Eq. ([3]) 571 
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where 11(0) comes from the subtraction to the corresponding dispersion relation due to the 
degree of divergence of the correlation function of the scalar glueball, and 
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where Gl° nt (so; s, r) is the contribution of continuum, 
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and Qj? (s, r) is defined as 
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with the Borel transformation B being defined by 
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The Gaussian sum rule emphasizes only the contribution of the hadron state considered, and 
suppresses the background exponentially (according to the Gaussian distribution). Recall 
that the Laplace sum rule stress only the contribution from the lowest state, and suppresses 
the other contributions exponentially (according to the exponential distribution). 
For k = — 1, 0, and 1, a straightforward but tedious manipulation leads to 



i 



\/Atxt Jo 



dss exp 



s — s) 



At 



X 



(oo - Tr 2 a 2 ) + 2ai In ( — — j + 3a 2 In 2 f — 



1 



+ 



+ 



V47TT Jo 

1 l"X} 



\/AtXT Jo 



dss exp 
(is exp 



s-s 



At 



s-s 



\fl 



TXT 



: exp 



V4 



: exp 



TXT 



s 
"47 

il 

~ At 
7 



4r 

C4tt 2 np~ 2 -— 

s 4r 

12 



167r 3 np 4 J 2 (pVs) F 2 (pv^) 
na 2 (C 2 p 2 ) 

(Coa s fnx + Ci« 2 n + 2 7 7r 2 n) 

5 \ TX I S 



-ln(v^p 2 ) + --]n2 



(25) 



^QCD 




16tt%p 4 J 2 ipy/s) Y 2 (pVs) 
na 2 (C 2 p 2 s - C 3 ) 



(26) 



S? CD ( S » 



V4VTT 7o 



<iss 3 exp 



.«-« 



X 



(a — vr 2 a 2 ) + 2oi In I — ^ I + 3a 2 In 2 I — ^ 



4r 
s 



v47rr Jo 
1 



y/4lTT Jo 



dss 3 exp 
dss exp 



s-s 



4r 



s-s 



At 



16vr 3 np 4 J 2 (pVi) ^2 (pVa) 
na 2 (C 2 p 2 s - C 3 ). 



V. NUMERICAL ANALYSIS 



(27) 



Now, we specify the input parameters in the numerical calculation. The color and flavor 
numbers are taken to be N c = 3 and Nf = 3, respectively. The expressions for two-loop 
quarkless (Nf = 0) running coupling constant a s (Q 2 ) at renormalization scale p 58l. |59| 
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and for the three-loop running coupling constant with three flavors (Nf = 3) 



(28) 
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are used, where the central value of the MS QCD scale A is taken to be 120 MeV, and 

L = lnfeV Li = ln 2 L-lnL-l, 

A, = J [11 - |iV>] , ft = £ [102 - f N f ] , 

a _ 1 T2857 5033 at- , 325 at- 1 
P2 - 43 L~2 IT-^/ + -54 ^f\ ■ 
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Recall here that research on the renormalization group improvement for Gaussian sum rules 



indicates that fi 2 = \fr 60j. The subtraction constant 11(0) is fixed by QCD low-energy 
theorem [3j 

n(0) = — (a s G 2 > ~ 0.6 GeV 4 . (31) 

9 

The values of the average instanton size and the overall instanton density are adopted from 



the instanton liquid model |42 



n = lf m - 4 = 0.0016 GeV 4 , p= -fm= 1.689 GeV" 1 . (32) 

Finally, the mass of the neutral pion is taken from the experimental data, i.e. m n = 135 
MeV. 

In order to measure the compatibility between both sides of the sum rules (1191) realized 
in our numerical simulation, we introduce a variation, 5, defined by 

1 " [L(rQ - R(n)] 2 

N^ \Hn)R(n)\ ' {66) 

where the interval [r m i n , r max ] is divided into 100 equal small intervals, N = 101, and L{ji) 
and Riji) are left-hand side and right-hand side of Eq.(TT9l evaluated at Tj. 

To determine the values of the resonance parameters in Eq. ( TT5l) . we match both sides 
of sum rules Eq. (fT9"|) optimally in the fiducial domain. In doing so, the value of s should 
approximately be set to be the corresponding mass squared of the resonance m 2 . To suppress 
the continuum contribution, we require s < m 2 . The conditions for determining the value of 
So are: first, it should be grater than m 2 ; and second, it should guarantee that there exists a 
sum rule window for our Gaussian sum rules. We note that the upper limit r max of the sum 
rule window is determined by requiring that the contribution from the continuum should be 
less than that of the resonance 

Qr'iso; s, r max ) < £ fe QCD (s ; s, r max ), (34) 

while the lower limit r m i n of the sum rule window is obtained by requiring the contribution of 
pure instantons to be greater than 50% of Q^ (sq; s, r), because such classical contributions 
should be dominant in the low-energy region. Moreover, to require that the multi-instanton 
corrections remain negligible, we simply adopt a rough estimate 



With these requirements, the figurations and numerical results are given below. 

For the case of quarkless QCD with the lowest resonance in the ++ channel, we adopt 
the isolated lowest resonance model (1151) for the spectral function, the optimal parameters 
governing the sum rules are listed in Table [H The corresponding curves for the left-hand 
side and right-hand side of (lT9"j) of k = — 1, 0, and +1 are displayed in Fig. [TJ where the solid 
lines are the right-hand side (QCD) of Eq. (!T9"j) . and the dashed lines are the left-hand side 
(HAD) of Eq. f|T9|) . and the dotted lines are for the right-hand side (QCD) excluding the 
contribution of interactions between the instantons and the quantum gluons (the same for 
hereafter). Taking the average, the values of the mass and width of the lowest ++ scalar 
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FIG. 1: The curves for the left-hand side and right-hand side of the Eq. (|19j) for quarkless QCD 
with only the lowest resonance considered. The solid line denotes the right-hand side (QCD), 
dashed line for left-hand side (HAD), and dotted line for the right-hand side (QCD) without the 
interaction contribution of the Gaussian sum rules (1191) . 
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TABLE I: The fitting values of the mass m and width T of the lowest ++ scalar glueball, and of the 
parameters A and f,so, [T m j n ,T max ], and 5 characterizing the couplings to the lowest resonance, 
the continuum threshold, the sum rule window, and the compatibility measure for finite-width 
Gaussian sum rules (|19p of k = — 1, 0, and 1 in quarkless QCD for a given s, the real part of the 
complex q 2 (q 2 = s + iQ 2 ). 

IjGcV 2 ) k m(GeV) r(GeV) A(GcV) /(GcV) s (GeV 2 ) [r min ,r max ](GeV 4 ) S 

1.49 2 ^1 L49 O04 1.481 1.495 4~40 [0.75-2.0] 2.40 x 10" 5 

1.50 2 1.50 0.15 1.523 1.534 4.65 [1.0-2.6] 7.79 x 10~ 6 

1.50 2 1 1.51 0.09 1.517 1.530 4.68 [0.8-1.5] 2.80 x 10~ 6 



glueball living in quarkless QCD, and the corresponding optical fit parameters are predicted 
to be 

m = 1.52 ± 0.18 GeV, T = 0.2 ± 0.15 GeV, 

/ = 1.47 ± 0.13 GeV, s = 3.8 ± 0.9 GeV 2 . (36) 

where the errors are estimated from the uncertainties of the spread between the individual 
sum rules, and by varying the phenomenological parameters, A and (a s G 2 ), appropriately 
away from their central values (the same for hereafter) 

A = 120-200MeV, (37) 

(a s G 2 ) = 0.6 -0.8 GeV 4 . (38) 

For quarkless QCD, there is only one well-defined scalar bound state of gluon suggested 
by lattice QCD and also from our investigation just described above. Including quarks 
enhances the difficulty of the task since many states possessing the same quantum numbers 
may be present in the correlator. Even so, at the first step as comparison, we still consider 
only the lowest resonance in the ++ channel for the case of QCD with three massless quarks 
as usual for the scalar glueball mass m of about 600 MeV, because that resonance / (600) 
may be considered to be well isolated. The optimal parameters governing the sum rules are 
listed in Table [TT1 The corresponding curves for the left-hand side and right-hand side of 
(TT9~|) of k — — 1, and, +1 are displayed in FigO Taking the average, the values of the mass 
and width of the (probable isolated) lowest ++ scalar glueball in the world of QCD with 
three massless quarks, and the corresponding optical fit parameters are predicted to be 

m = 1.54 ±0.17 GeV, T = 0.23 ± 0.13 GeV, 

/ = 1.64 ± 0.14 GeV, s = 3.8 ± 0.9 GeV 2 . (39) 

The above one isolated lowest resonance assumption is questioned from the admixture 
with quarkonium states, and from the experimental data that three ++ scalar states are 
around the mass scale of 1500 MeV [namely / (1370), / (1500), and / (1710)]. Therefore, we 
adopt the three-resonance model for the phenomenological side of the sum rules for the case 
of QCD with three massless quarks in the ++ -channel when m is above 1 GeV; the optimal 
parameters governing the sum rules are listed in Table IHIl The corresponding curves for 
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FIG. 2: The curves for the left-hand side and right-hand side of the Eq. (]19p for three-flavor QCD 
in chiral limit with only the lowest resonance considered. The solid line denotes the right-hand 
side (QCD), dashed line for left-hand side (HAD), and dotted line for the right-hand side (QCD) 
without the interaction contribution of the Gaussian sum rules (1191) . 
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FIG. 3: The curves for the left-hand side and right-hand side of the Eq. (]19p for three-flavor QCD 
in chiral limit with the lowest three resonances considered. The solid line denotes the right-hand 
side (QCD), dashed line for left-hand side (HAD), and dotted line for the right-hand side (QCD) 
without the interaction contribution of the Gaussian sum rules (1191) . 
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TABLE II: The fitting values of the mass m and width V of the lowest ++ scalar glueball, and 
of the parameters A and /, sq, [r m i n ,r max ], and 5 characterizing the couplings to the lowest single 
resonance, the continuum threshold, the sum rule window, and the compatibility measure for finite- 
width Gaussian sum rules (|19p of k = — 1, 0, and 1 in QCD with three massless quarks for a given 
s, the real part of the complex q 2 (q 2 = s + iQ 2 ). 



s(GeV 2 ) 


k 


m(GeV) r(GeV) A(GeV) /(GeV) s (GcV 2 ) [r min ,r max ](GeV 4 ) 6 


1.45 2 


-1 


1.52 0.10 1.590 1.602 4.5 [0.5-2.0] 6.87 x 10" 5 


1.50 2 





1.53 0.10 1.635 1.646 4.7 [1.5-5.0] 1.34 x 10~ 5 


1.52 2 


1 


1.53 0.10 1.663 1.674 4.7 [1.5-2.7] 3.04 x 10" 5 



TABLE III: The fitting values of the mass m and width T of the lowest ++ scalar glueball, and 
of the parameters A and /, so, [r m i n , r max ], and 5 characterizing the couplings to the three closely 
located resonances, the continuum threshold, the sum rule window, and the compatibility measure 
for finite- width Gaussian sum rules (|19p of k = — 1, 0, and 1 in QCD with three massless quarks 
for a given J, the real part of the complex q 2 (q 2 = s + iQ 2 )- 



s(GeV 2 ) 


k 


m (GeV) r(GcV) A(GeV) /(GeV) 


s (GcV 2 ) [r mil 


„r max ](GeV 4 ) 


S 






1.37 


0.30 


0.950 


0.983 








1.35 2 


-1 


1.50 
1.71 


0.10 
0.14 


1.510 
1.100 


1.523 
1.125 


4.2 


[2.0,8.0] 


2.90 x 10~ 4 






1.37 


0.30 


0.950 


0.983 








1.35 2 





1.50 
1.71 


0.10 
0.14 


1.583 
1.100 


1.595 
1.125 


4.0 


[3.0,8.0] 


6.42 x 10~ 6 






1.37 


0.30 


1.100 


1.125 








1.35 2 


1 


1.50 
1.71 


0.10 
0.14 


1.600 
1.100 


1.612 
1.125 


4.0 


[3.0,8.8] 


5.13 x 10~ 6 



the left-hand side and right-hand side of ( TT91) of k = — 1, 0, and +1 are displayed in Fig. 
Taking the average, the values of the widths of the three lowest ++ scalar resonances in the 
world of QCD with three massless quarks, and the corresponding optical fit parameters are 
predicted to be 

m = 1.37 ± 0.06 GeV, T = 0.30 ± 0.10 GeV, / = 1.10 ± 0.13 GeV. (40) 

for the resonance / (1370), and 

m = 1.50 ±0.10 GeV, T = 0.10 ± 0.06 GeV/ = 1.60 ± 0.11 GeV. (41) 

for / (1500), and 

m = 1.71 ±0.11 GeV, T = 0.14 ± 0.08 GeV/ = 1.10 ± 0.14 GeV. (42) 
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for / (1710). 

Figures{T]l3] show the consistent match between the both sides of Eq. (TT9l) for k = — 1, 0, 
and 1, respectively, with the fitting parameters. The matching between both sides of the 
sum rules is very well over the whole fiducial region with a very little departure. 

These results are in good accordance with the experimental data of / (1500), m = 1505±5 
MeV, T = 109 ± 7 MeV Q, and the sum rule calculation of Ref. 0, m = 1.53 ± 0.2 GeV, 
/ = 1.01 ± 0.25 GeV (see Table HV|) . We do not calculate the higher moments sum rule, 
because for k > 1, the continuum contributions become very large. 



TABLE IV: The numerical results from other methods. 



Methods 


Mass (GeV) 


Width 


(GeV) 


Coupling (GeV) 


■so (GeV 2 ) 


leferences 


GSR 


0.8- 1.6 


0.4 - 0.6 




2.3 


2JJ 


QLQCD 


1.3- 1.7 


0.035 - 


0.873 






li-18| 


QSSR 


1.4 ±0.2 






2.56-2.61 




g|| 


LSR 


1.25 ±0.2 
1.53 ±0.2 
1.52 ±0.2 






1.05 ±0.1 
1.01 ±0.25 
0.39 ±0.145 


5.0 ±0.1 
5.0 ±0.1 
4.2 ±0.2 




model 


1.42- 1.5 

1.666 

1.633 
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VI. DISCUSSION AND CONCLUSIONS 

The properties of the ++ scalar glueball are examined in a family of the finite-width 
Gaussian sum rules. The correlation function is calculated in a semiclassical expansion, a 
well-defined process justified in the path-integral quantization formalism, of QCD in the 
instanton background, namely the instanton liquid model of the QCD vacuum. Besides 
the contributions from pure gluons and instantons separately, the one arising from the 
interactions between the classical instanton fields and the quantum gluon ones are taken 
into account as well. Instead of using the usual zero-width approximation for the spectral 
function of the considered resonances, the Breit-Wigner form for the resonances with a 
correct threshold behavior is adopted. With the QCD standard input parameters, three 
Gaussian sum rules with the k = — 1, 0, and 1-th moments are carefully studied. 

For the quarkless QCD, we have, in fact, changed the value of s, and found that the value 
of the mass of the lowest resonance is approximately proportional to s, and the value of sq 
arrives almost at its maximum for s lying between 1.50 GeV 2 and 1.70 GeV 2 , where the 
couplings to the state are almost the same for different k, and the corresponding widths 
become small and stable. We have only shown the situation with the optimal compatibility. 
We note here that the value of 5 for m = 650 MeV is one or two-orders lower than the 
optimal one, and the values of / are not coincident for different k, so that the mass scale 
of the lowest ++ scalar glueball may not be lower than 1 GeV. The mass and width of the 
lowest glueball without quark loop corrections are predicted in Eq. (1361 ). 
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For QCD with three massless flavors and by considering only single scalar resonance, the 
same behavior with respect to the changing of s appears. Namely, the value of the mass of 
the resonance is approximately proportional to s, and the value of so arrives almost at its 
maximum for s lying between 1.35 GeV 2 and 1.70 GeV 2 , where the couplings to the state 
are almost the same for different k, and the corresponding widths become small and stable. 
When s = 1.70 2 GeV 2 (the situation with optimal compatibility), all physical parameters are 
almost the same for different k. We note here that the value of S for m = 650 MeV is nearly 
one order lower than the optimal one, and the values of / are not coincident for different k 
as in the case of pure QCD, so that the mass scale of the lowest ++ scalar glueball, even 
in the world of QCD with massless quarks, may not be lower than 1 GeV. The mass and 
width of the lowest glueball with quark loop corrections under the assumption of no mixture 
between glueball and qq state are predicted in Eq. fl39|) . 

For QCD with three massless flavors and by considering three closely located ++ scalar 
resonances / (1370), / (1500), and / (1710) (namely both s and the resonance masses are 
given input parameters), the behavior with respect to the changing of s is changed. Namely, 
the values of the widths and Sq for the three resonances remain as almost invariant, and 
only their couplings slightly increase when increasing s. The optimal compatibility is arrived 
at s — 1.35 2 GeV 2 , as shown in Table HTT1 and Figj3j The widths of the lowest three ++ 
resonances coupled to the glueball current O s are predicted in (jlOj) , from which we can read 
off the corresponding couplings f 3 to the three resonances /o(1370), /o(1500), and /o(1710) 
with masses 1.35, 1.47, and 1.70 GeV as 

0.95 ±0.47 GeV 3 , 3.92 ± 0.85 GeV 3 , 1.42 ± 0.51 GeV 3 , (43) 

respectively. 

In summary, we may conclude that, first, any ++ scalar resonance below 1 GeV, such 
as /o(600), contains almost no component of the scalar glueball; second, the values of the 
mass and decay width of the ++ resonance, in which the fraction of the scalar glueball 
state is dominant, are m = 1.50 ± 0.10 GeV and T = 0.10 ± 0.06 GeV, respectively, and 
the value of its coupling to the corresponding current is / 3 >m 2 = 3.92 ± 0.85 GeV 3 ; third, 
the fractions of the scalar glueball contained in the other nearby ++ scalar resonances, 
/o(1370) and ,fn(1710), are also appreciable. These are not only compatible with lattice 
QCD simulation [14Ml8| and other estimations 2314251. 129. 6lL but also in good accordance 



with the experimental data of the low scalar resonances 34-136 

It is also remarkable that the three Gaussian sum rules lead to almost the same results, 
a consistency between the subtracted and unsubtracted sum rules is very well justified. 
We note that we have not been working within the mixed scheme, namely with including 
condensates, and in the same time, adopting the so-called direct-instanton approximation, 
but simply with a self-consistent framework, a quantum theory in a classical background, 
without the problem of double counting. In this aspect, our results further justified the 
instanton liquid model for QCD among many other justifications. 

In our semiclassical expansion, the leading contribution to the sum rules comes from 
instantons themselves, especially in the region below the threshold sq. It is the amount of 
this contribution that determines the low bound of the sum rule window. This means that 
the nonlinear configurations of gluons have a dominant role with respect to the quantum 
fluctuations in the low-energy region. 

The contribution of the interactions between the classical instanton fields and 
gluon ones, considered in this paper but neglected in earlier sum rule calculations 

16 



l q uaniuE 
2a|2l|22 



29|, is in fact not negligible. To the contrary, its amount is approximately double or even 
triple that from the pure quantum fluctuations in the whole fiducial domain, expected from 
a view point of the semiclassical expansion. Moreover, it is obviously seen from Figs. |T](3] 
that, without taking the contribution from the interactions between instantons and quan- 
tum gluons into account, the departures between Gk ad (s ;s,T) and Qjf (s ;s,r) without 
interaction become large, and all three Gaussian sum rules become less stable, and thus less 
reliable. 

Finally, it should be noticed that the imaginary part of the instanton contribution is an 
oscillating, amplifying and nonpositive defined function, and so is the imaginary part of the 
correlation function. This property which is a fatal problem for the QCD sum rule calculation 
with the instanton background, may make the contribution of continuum too large to be 
under control. Hilmar Forkel introduced a Gaussian distribution for the instanton to get rid 
off this trouble, and obtained a smaller ++ mass scale: 1.25 ± 0.2 GeV 23j compared to 
the earlier result 1.53 ±0.2 GeV [24J . We did not use this Gaussian distribution, but simply 



chose a smaller fitting parameter to avoid this problem. 
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